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Abstract 

4^ Wc consider the stellar production of vector states V within the minimal model of "dark 

photons" . We show that when the Stiickelberg mass of the dark vector becomes smaller than 
Cn| plasma frequency, the emission rate is dominated by the production of the longitudinal modes 

^ of V, and scales as K^my, where k and my are the mixing angle with the photon and the 

OO mass of the dark state. This is in contrast with earlier erroneous claims in the literature that 

^ the emission rate decouples as the forth power of the mass. We derive ensuing constraints on 

^ the (k, my) parameter space by calculating the cooling rates for the Sun and red giant stars. 

O We find that stellar bounds for my < 10 eV are significantly strengthened, to the extent that 

^ all current "light-shining-through- wall" experiments find themselves within deeply excluded 

regions. 

>< 



1 Introduction 



The Standard Model of particles and fields (SM) can be naturally extended by relatively light 
neutral states. Almost all possible ways of connecting such states to the SM have been ex- 
plored, and several of such ways stand out as the most economical/natural. One of the most 
attractive possibilities is the so-called "hypercharge portal", or "kinetic mixing" portal that 
at low energy connects the electromagnetic current with another massive photon-like state 
[1]. This model has been under intense scrutiny in the last few years, both experimentally 
and observationally. The interest to this model is fueled by attractive (yet speculative) pos- 
sibilities: the dark vector can be a promising mediator of the dark matter-SM interaction [2j , 
or form super- weakly interacting dark matter itself [21 Dark vectors were proposed as a 
possible solution to the muon g — 2 discrepancy [5] , and have been searched for (so far with 
negative results), both at high energy and in medium energy high- intensity particle physics 
experiments. 

The region of small vector masses, my < eV, can also be very interesting. On the 
theoretical side, there are speculations of dark photons contributing to dark matter (via an 
initial condensate-like state) [B] and dark radiation [7j. But perhaps more importantly, there 
are some hopes for the terrestrial detection of dark photons. So, far several avenues have 
been proposed: one can attempt observing a " visible-dark- visible" oscillation chain in "light- 
shining-through- wall" experiments (LSW) |8j. The quanta of dark photons emitted from the 
Sun can be searched for with "helioscopes" [9], neutrino- pjij and dark matter experiments 
p!T| [12]. Some of these exciting possibilities have been summarized in the recent review [13] . 
We will refer to all proposals and experiments aimed at detection of dark vectors, produced 
astrophysically or in the laboratory, as direct searches. 

At the same time, it is well-known that for many light (my < keV) and weakly-coupled 
exotic particles the astrophysical constraints are often far stronger than direct laboratory 
constraints The astrophysical constraints are very important for the dark vectors as 
well, as they determine a surviving fraction of the parameter space that can be explored in 
direct searches. The most important limits to recon with are the constraints on the emission 
of dark vectors from solar luminosity, from the red giants, neutron star and supernovae 
cooling rates. 

To date, the only in-depth analysis of astrophysical bounds on sub-keV dark vectors 
was performed by Redondo in [9]. In this paper we re-assess these bounds, provide correct 
calculations for the dark photon emission rates, and strengthen the astrophysical bounds in 
the LSW region by as much as ten orders of magnitude. Our findings significantly reduce 
the parameter space available for the direct searches and affect or completely change the 
conclusions of many papers written on this subject. In a separate forthcoming publication 
we will address new limits imposed by the most advanced WIMP detectors on the solar 
emission of dark vectors [T5] . 

This paper is organized as follows. The next section introduces the minimal model of 
the dark vector, and explains the main scaling of its production rate with my. Section 3 
contains technicalities of the in-medium production of the dark vector. Section 4 contains 
practical formulae for the stellar emission rates, in application to the Sun and red giants. 



1 



and sets the constraints on the mass-mixing parameter space. We reach our conclusions in 
section 5. 



2 Dark photon production, in vacuum and in a medium 

The minimal model of "dark vectors" extends the SM gauge group SU(3)c x SU{2)l x U(l)y 
by an Abelian factor U(l)y. Kinetic mixing of the hypercharge field strength with 
the field strength V^^, of U(l)y links the SM to the new physics sector, while SM fields are 
assumed to be neutral under U(l)y. We are interested in processes far below the electroweak 
energy scale, for which the relevant low-energy Lagrangian takes the form 



Here F^^ = d^A^ — d^A^ is the photon field strength and is the "hidden photon" (also 
known as "dark vector", "secluded vector", "dark photon" etc — an equivalent set of names). 
The coupling of A^ and V"^ is regulated by the kinetic mixing parameter k, redefined in an 
appropriate way to absorb the dependence on weak mixing angle. For all calculations in this 
paper we use k -C 1, and consider only leading order terms in the mixing angle. Finally Jl^^ 
is the usual electromagnetic current with electric charge e < 0. 

It is important to comment on the origin of my in ([T|. The simplest possibility is that my 
is a Stiickelberg-type mass. Because of the conservation of the Abelian vector current, my 
remains protected against sensitivity to UV scales, and such a model is technically natural 
even with very small my. An alternative generic possibility is a new scalar field (s) charged 
under U(l)y that develops a vacuum expectation value that Higgses the hidden group. This 
introduces a new interaction term of the physical hidden Higgs with vectors, g'myh'V^, as 
well as h' self-interaction (see e.g. [IS]). It is well understood that in the limit of my and 
mhi small compared to all energy scales in the problem, the production of dark sector states 
is dominated by the dark Higgsstrahlung [31 [TH|, or equivalently, by the pair-production of 
the U(l)y-charged Higgs scalar fields. Importantly, this process is insensitive to the actual 
mass of my in the small mass limit, and schematically 



where we show only the dependence on dark sector parameters, leaving the SM part of the 
V + h' production process completely general; a' = {g'Y /{An) is the square of the coupling of 
dark Higgs to V^. For sub-keV dark vectors and Higgses, all previously derived constraints 
on " millicharged particles" apply [T3|, and limit the Kg' combination to be below ~ 10^^'^. 
The technical reason for not having any small my suppression of the rate ^ despite the 
interaction term g'myh'V^ being proportional to my is of course tied to the production of 
the longitudinal modes oiV mV + h' final state. 

The models with the hard {i.e. Stiickelberg) mass my behave differently as the production 
rate of dark vectors has to decouple in the small my limit. The easiest way to see this is to 
restrict the interaction terms in ([T]) to on-shell V^, using d^V^ = and to leading order in 



£ = --F^ - -V^ --F + —V + eJ^" A 



(1) 



RatesM^v+h' oc a'/t^(my)°. 



(2) 
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Figure 1: Illustration of the dark photon emission process by the electromagnetic current. 
n, ^^V^"' = -m\V\ so that 

Ant = -^FM-^'^" + e^c^^M L,^, = -Km\A^V^^er,^A^. (3) 

This expression is of course explicitly gauge invariant under — > + d^x due to the 
current conservation and on-shellness of conditions: 

a^Ji^, = 0; 9^r^ = 0. (4) 

The appearance of in the coupling of and shows that two sectors are decoupled 
in my = limit. The most important question in considering the production of states 
is the scaling of the production rate with my, in vacuum and inside a medium. The exist- 
ing literature on the subject |9] and its subsequent follow-up papers claim that in-medium 
production decouples as Rate su^v n'^rriy in the small mv limit. This inference is wrong. 

To demonstrate our point we consider a generic production process i ^ f + V due to 
(|3|, where i, f are any initial, final states of the SM particles. A schematic drawing of 
such a process is shown in Fig. [TJ Without loss of generality we assume that V is emitted 
in 2;-direction, so that its four-momentum is given by (a;,0,0, with uJ^ — k'^ = my. 
Moreover, we assume that the energy of the emitted V is much larger than its rest mass, u 3> 
mv- Three polarization states can be emitted: two transverse states Vt with polarization 
vectors = (0,1,0,0,) and (0,0,1,0), and one longitudinal mode with polarization 
vector e-^ = my^{\k\,Q,Q,uj). In all cases = —1 and e^k^ = 0. 

We include a boundary-free medium via some conducting plasma, characterized by the 
plasma frequency Up. We consider two regimes, [almost] vacuum: Up -C my <C oo, and 
in-medium: my Up u. The choice of \k\,u ^ Up is not essential, and we consider all 
ranges of u in the next section. The matrix element for the production process induced by 
(|3| is given by 

Mi^f+v^^^^ = Km\ \eJ,r.^^n (v4^ N') e^^^), (5) 

where (A^, A'^) stands for the photon propagator with input momentum /c^, and [eJ^j^]/* is 
the matrix element of the electromagnetic current. We disregard various my-independent 
phase factors and normalizations, as our goal in this section is to only consistently follow 
the powers of my. 
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For convenience, we fix the photon gauge to be Coulomb, VjAj = 0, but same results 
can be obtained in other gauge, of course. The photon propagator for the production of the 
transverse modes is given by (see, e.g. [18]): 



(A,, A^) = $ = -> 5^ X ( >> (6) 



UJ 



2 _\^2 m\, - uj^ ^ y -Up at my < Up 



where Sji is the projector onto transverse modes. Because of the existence of two different 
regimes for the transverse modes of the photon, the amplitude of the V'^ production is 
my-independent in vacuum, and m^-suppressed in medium. 

The production of the longitudinal modes in Coulomb gauge is mediated by the {Aq, Aq) 
part of the propagator, li k u ^ Up, this propagator is unaffected by the medium [18], 



(Ao,Ao) = ^. (7) 

The difference in the behavior of the two parts of the photon propagator can be readily 
understood. While the propagating transverse modes of the photon are close to mass shell, 
and plasma effects can affect them easily no matter how large u is, the non-propagating 
piece (^40,^0) cannot be changed in a causal manner if \k\ is larger than e.g. the inverse 
distance between particles in the plasma. Plugging the expression for the propagator 
together with the explicit form of into (p|), we arrive at the conclusion that the amplitude 
for the production of the longitudinal mode is proportional to the first power of my. 

Collecting all relevant factors in one expression, we get the all-important my scalings for 
the production of both transverse and longitudinal modes: 

Rate oc I '^^ vacuum, my > cUp, 

sm^Vt K^niyUp'^ in medium, my ^ Up. 

Performing the same estimate for the production of the longitudinal mode, (and keeping 
UJ ^ Up) we get 

Rate sm^Vl 0*^ K^myco^'^, both in vacuum and in medium. (9) 

Notice the quadratic, and not quartic dependence on the dark photon mass in ([9]), in contrast 
with conclusions of Ref . [9] . The my-scaling of the production rate for the transverse modes, 
Eq. ([s]), is of course a well-known result in the "dark photon" literature. It exhibits rriy 
decoupling at small my, but as it turns out, this was incorrectly extended to the production 
of the longitudinal modes. (Thus we also learn that "nature does not like to skip an order," 
and compared to the Higgsed case of an 0(m^)-rate ^ the emission of Stiickelberg vectors 
occurs in lowest possible 0{my) order.) We also note in passing that for some processes, 
only the production of of the longitudinal modes of is actually possible. A well-studied 
process is the — )■ tt^V decay [3], which is fobidden for the transverse modes, but has the 
expected my scaling at small my for the longitudinal modes. 

The correct scaling ^ will bring about momentous change in all estimates of stellar 
cooling rates. For the favorite LSW region with my ~ 10"'^ eV, one should expect that the 
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previous literature underestimates the solar cooling rates by as much as my/Up ~ 10~^°, 
and correspondingly, one expects the tightening of the constraints by the same large factor, 
once the production of the longitudinal mode of the dark vector is treated correctly. In the 
two subsequent sections we perform this analysis in some detail. 



3 Plasma production rate of dark photons 

Inside a medium, the propagation of the electromagnetic field is determined by the electro- 
magnetic polarization tensor 11'^'^ = e^(J^, Jem)- Due to k^Jj^^ = where k = {u, k) is the 
four-momentum flow inside the polarization tensor, 11'^'^ can be parameterized as 

n'^^ = Ht 5^ e^fe^r + n^e^'e^^ , (10) 

j=l,2 

where e^'*^ are the transverse and longitudinal polarization vectors of a vector boson with 
momentum k. In particular, 

6^= , ^ -Mkluji^) . ill) 



\k\ 



In the Coulomb gauge, the propagator of the electromagnetic field can be written as 

1 / y^* y 

iA\A^) = ^ U^^-^ 

Lo^ -\k\'' -UtK k^ 

(^^^°) = . (12) 



where eg = \k\/ \Jk^ki^ has been used. Notice that while the definition of IIt- is usually 
uniform across the literature, the definition of fl^, varies, and e.g. in Ref. [18] it is defined 
differently, II^^'^^' ^ = j^L^^ n^**^"^ work^ With the explicit expression of the photon propa- 
gator, the matrix element for the dark photon emission in Eq. ([s]) can further be written 



as 



Km 



2 



Illy — li-T 



M.^nv. = J^'^^AeJ^h^ ■ (13) 



V 



^L\k\ 



Using the condition k^^J^^ = and Eq. (11), it is easy to show that 



■^em^O ~ _2 '^ern^u ■ (1^) 



\k 



2 



m 



2 era fi 
V 
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Therefore, Eq. (13) can be further simphfied to 



2 

M..^f^vr,. = --2^^^^[eJe^]/.eJ^ . (15) 

niy — iiT,L 

The general expression for Re IIt and Re 11^ can be found in Ref . . Thus we find that the 
emission of dark vectors is given by the vacuum matrix element for the emission of massive 
photons, [eJ^]/ie^'^, with fiducial photon mass = my and multiplied by the effective 
mixing angles, defined according to, 

2 4 

2 K, rUy . , 

^T,L = 7—2 P„TT A. , /T^TT ■ (16) 



m 



V 



At finite temperature T 7^ 0, the imaginary parts of IIt.l are related to the rates at which 
the respective distribution functions approach equilibrium fT9]. Detailed balance equation 
allows this to be expressed exclusively in terms of the absorption rate T^^^^^ of in-medium 
transverse (longitudinal) massive photons, 

ImUrA^, \k\) = -u{l- e-^l^) \k\) . (17) 

It will be convenient to express T^^^^-j in terms of a differential production rate dr^°j^/ [dudV) 
per frequency interval and volume. To this end we write a Boltzmann-type equation for the 
photon distribution function dnT,L/du, 

driT^L ^ d'^rf 1 dnT,L /-.on 

diodt dudVl-e-^/T du ^ ' 

Using detailed balancing and noting that the equilibrium photon distribution is given by 
driT^L/doj = 5'T,Li^|^|/(27r^)/(e^'^/"^ — 1), where gr = 2 and gL = 1, are the degeneracies for 
the transverse and longitudinal modes respectively, this yields, 

pabs _ 2^ T,L ^IT /.q^ 



k\ dudV 



Therefore, according to Eqs. (15), (16), we conclude that the differential production rate 
of transverse and longitudinal modes of dark vectors, dT^°^'^ / (dudV) , can be written as 

^pprody ^ ^ ^pprod 

Current conservation demands that both RcIIl and Imlli are proportional to my. 



Therefore, k| will remain my-independent in the small my limit as can be seen from (16) 
when expressed as 



'^i ~ / N 2 / N 2 ' C^^) 

RcHr \ I / ImHr 



1-^ + 
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From Eqs. (16) and (21) it is obvious that in the region my ^ RellT, scales as rriy, 
whereas k| scales as my. However, r^™*^ scales as my due to the current conservation. 
Therefore, the production rate of the transverse modes is suppressed by my in the small 
mv limit whereas the production rate of the longitudinal mode is only suppressed by my, 
in agreement with our qualitative discussion in the previous section. 

Finally, for a finite size thermal system, the radiation power into the dark photons can 
be written as 

P = [dV [dco (4^^ + .1"^^] . (22) 
./ / I ^ dudV ^ dudV / ^ ^ 



In the next section, we use (22) to calculate the V^-emission rate as a function of my. 



4 Solar luminosity and red giant stars constraints on dark vectors 

Inside the sun, the electrons are non-relativistic and non-degenerate. In this limit, the leading 
order terms in a T/me expansion of RcIIt^l can be written as 

ReHr = , ReUL = ^1 - j ^ (23) 

and the plasma frequency is given by 

^l = — = , 24 

^ me me 

where rie is the number density of electrons. 

The production rate of the transverse modes is thoroughly calculated in Ref. [9], which 
we have checked and agree with. As already mentioned multiple times, the production of 
longitudinal modes is treated incorrectly in [U]. We believe that one can trace the error 
to n^'^^''^ defined as — in Eq. (8) of Ref. [H]. This is unphysical, as 11^ generated 
by the plasma must vanish in the ng — )■ limit. We encounter the same problem in the 
corresponding chapters of Ref. p^. We can only speculate that this error crept in from 
an inaccurate translation of the results from Ref. [18] obtained in the Coulomb gauge to a 
generic co variant gauge. 

The production rate of the longitudinal modes has contributions both from bremsstrahlung 
and from Compton scattering. Inside the sun and a red giant star the bremsstrahlung pro- 
cesses are dominant. In the non-relativistic and non-degenerate limit the cross section for 
the bremsstrahlung process to produce the longitudinal massive photons at fc^ = my can be 
written as 

do -y^Z'^f^mi^^ uj'^ — my ^ v 

= 3{2nrmlu'v^ y - _ 2u;/m, 

where v is the velocity of the incoming electron. Then, appropriately averaging over the 
electron energy distribution, the production rate can be written as 



(25) 



dTl'°'' 



dVdu 



brem " ^ 



3m2 a;4 V \ nT ' \ y T 
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Figure 2: Upper limits on the kinetic mixing parameter k, vs. the mass of the dark photon my. 
The sohd blue and red curves are the total constraints on the dark photon parameter space 
from the sun (blue) and red giants (red). The dashed curves show constraints from retaining 
only the longitudinal resonance contributions. For comparison, the dot-dashed curve shows 
the upper limit by considering only the contribution from transverse modes. The current 
bounds from the latest result of the LSW experiment by the ALPS collaboration is 
shown in gray, accompanied by the potential reach (yellow region) of the next generation of 
LSW experiments [26] . 
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where n^. is the number density of ions of charge — ZjC, and 



/■^ _ 2 

/(a) = / dx xe ^ log 

J a 



X + y/x"^ — a'^ 



X 



- ^Jx^ 



(27) 



In principle, inside the plasma, Eq. (25) should be modified to take into account the Debye 
screening effect characterized by the Debye screening length A^, = T/ie^rif.). However, 
inside the sun, a numerical study shows that the square of the typical momentum transfer of 
electrons ~ rn^T is much larger than 1/A|), and therefore this effect is of little importance 
for our level of rigour. 

Substituting the explicit dependence of RellL into k|. 



,2 



/ Im^L ^ 



2 



(2^ 



we observe that the production rate of the longitudinal modes of dark vectors reaches a 
resonance at w = ujp. Therefore, for my <^ Up the vector can be produced on resonance at 
any temperature. The resonant contribution to the radiation power per unit volume and per 
unit frequency can be written as 



dPr 1 K^mirOj: 



dVdu Att e'^p/'^ - 1 



x5{uj-Up) , (29) 



which is independent of the details of the production processes. Therefore, the resonant 
contribution to power emitted into longitudinal dark vectors can be written as 

PeLs-^^'ral^l dr ^^j)^' _ ^ . (30) 

where Rq is the radius of the sun. The corresponding results in Ref. [9] contain an additional 
factor of 2/7r X my/ojp under the integral, which we argue is wrong. 

We use the standard solar model BP05(OP) [20j to calculate the total power radiated 
into dark photons. Consistency of this model with observations, requires this additional 
" dark radiating power" be smaller than the actual measured solar luminosity [211 122] , which 
is Lq = 3.83 X 10^^ Watt. This sets the limit on the parameter space of the model, that 
we plot in Fig. [2j The solid blue curve shows the resulting constraint, and everything above 
this curve is excluded. For comparison, we also show the breakdown of the constraint by 
different contributions. The dashed blue curve shows the constraint resulting from retaining 
only the resonant contribution in the longitudinal modes. We can see that in the small my 
region (my < 0.1 eV), the dark radiation is indeed dominated by the longitudinal resonance 
contribution where the constraint can be simplified as 

«: X ^ < 1.4 X 10-" . (31) 
eV 

The dot-dashed brown curve shows the constraint by considering only the contribution 
from the transverse mode, which coincides with the solid blue curve in the large my region. 
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my > 10 eV, where we are also in agreement with [H]. However, in the small my region, 
the emission of T-modes gives only a subdominant contribution, and indeed the solid blue 
curve positions itself much below the brown dot-dashed curve. Therefore, our paper greatly 
improves the bounds on dark photons below the mass range of a few eV. 

For the red giant stars, we set the constraint by requiring that energy loss into dark 
photons should not exceed the nuclear energy generation rate. Thus, the radiation power 
into dark photons in the nuclear reaction region of the red giant should not exceed 10~^ 
Watt gram~^ 123] • The average temperature and density of the nuclear reaction region 
of red giant can be estimated as 

= lO^K , = 10^ gram/cm^ . (32) 

In this region, the electron gas can still be viewed as non-relativistic and non-degenerate Pl] . 
The solid red curve in Fig. [2] shows the upper limit on k by requiring the dark radiation 
power is smaller than the nuclear reaction power. The dashed red curve shows the upper 
limit by considering only the resonant longitudinal contribution. One can see that in the 
region my < 100 eV, the radiation is dominated by the longitudinal mode. The spike at 
my ~ 2 keV is caused by the resonant production of the transverse modes which appears at 
my = Up. Overall, we find that the solar luminosity provides a somewhat stronger constraint 
for all masses below several 100 eV. 

5 Summary and Discussions 

We have shown that the production rates of massive dark vectors V with the Stiickelberg 
mass my, coupled to the SM via the kinetic mixing portal, scales as rriy at small my due 
to the emission of the longitudinal modes of V. This drastically change the strength of the 
stellar constraints in the small my region. Thus, for the first time, and despite the large 
abundance of literature on dark photons, our paper sets correct stellar constraints on the 
dark photon parameter space in the whole mass range below a few eV. This turns out to 
be a region of special interest for LSW experiments. In Fig. 2, we show that even the most 
advanced ones (ALPS) [25] find themselves inside a deeply excluded region. Recalling that 
the signal of dark photons in LSW experiments scales two order of magnitude gap in k, 

at my ~ 0.01 eV between stellar constraints and LSW region translates into a required eight 
orders of magnitude improvement in sensitivity before LSW experiments become competitive 
with stellar bounds. Indeed, a large part of the sensitivity region for the next generation of 
LSW experiments, deemed reachable in Ref. [26], is also excluded. 

We conclude with several additional remarks pertaining to light dark photons: 

• The correct scaling with my, the m^-behaviour, is important not only for the produc- 
tion, but also for the detection of dark photons in the laboratory environment. Taking 
an example of dark matter detectors, made of some material with refractive index n, 
one can expect that in the region rriy <^ a;^|l — n|, the absorption of the transverse 
mode scales as my, whereas the absorption of the longitudinal mode scales as rriy. 
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Therefore, in this region one expect that the longitudinal mode will dominate the ob- 
served signal. The details of setting constraints on the solar dark photons with the use 
of the most sensitive low-energy threshold dark matter detectors will be addressed in 
our next paper [15j. 

• It has to be emphasized that in this paper we assume that, once produced, the dark 
photon freely escapes the stellar interior. It is possible that for large values of mixing 
angles, the energy loss process is quenched because of absorption. Determining whether 
such "islands" actually exist deep inside the excluded region, and if so whether they 
survive other constraints goes outside the scope of the present paper. 

• Stellar constraints derived in this paper have implications for the "dark CMB" [7j. We 
do not expect that the emission of the longitudinal modes in processes like 7r -|- e — )■ 
Vi+ein the primordial plasma will drastically change the estimates for the total energy 
density locked in V modes, because of the very small value for rie/n^ ~ 6 x lO"^*^. The 
combination of our new stellar constraints derived in this paper, and the constraints 
from spectral distortions of the normal CMB in the very small my region prevents 
generating large deviations of the effective number of relativistic degrees of freedom 
N^ff from its standard model value. 

• Finally, we can go back to the picture of the "Higgsed" dark photon production, Eq. 
(|2|, and ask the question of whether one can avoid strong constraints on k by choosing 
very small a', and by having the solar interior restore the symmetry in 11(1)^ sector, so 
that my = inside the sun. The symmetry restoration may occur for KLOp > my, which 
can have an overlap with the LSW regions of interest. Such an unusual "chameleonic" 
scenario deserves a special analysis. 
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